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Abstract 

By using suitable Ward identities, we show that it is possible to compute K°-K° 
mixing without subtracting the terms generated by explicit chiral symmetry break- 
ing present in Wilson-like lattice actions. The accuracy in the determination of the 
amplitudes is of 0(a), which is the best one attainable in the absence of improvement. 



PACS: 13.75Lb, 11.15.Ha, 12.38.Gc. 



1 Introduction 



A key ingredient in the study of CP violation in the Standard Model is the theoretical 
prediction of the K°-K° mixing amplitude. This involves the computation of the 
AS = 2 matrix element 

(K°\O AS = 2 \K ) ee \f\m\B K (1) 

of the operator 

Q AS=2 = Qi = _ 75) ^ )( -B^ (1 _ 7s)dB) ; (2) 

where s and d stand for strange and down quarks and A, B are colour indices. 

In general, important information on the physics beyond the Standard Model, such 
as various SUSY extensions (MSSM, NMSSM,...), can be obtained by studying AF = 2 
transitions (see |l]] and references therein for a discussion). For neutral kaons, such pro- 
cesses require, besides 0%, also the knowledge of the matrix elements of the operators 
(we adopt here the notation of ref . M ) 

2 = (s A (l- l5 )d A )(s B (l- l5 )d B ) 

3 = (s A (l- l5 )d B )(s B (l- l5 )d A ) 

4 = (s A (l - 7 5 )d A )(s B (l +K)d B ) (3) 

5 = (s A (l- l5 )d B )(s B (l + l5 )d A ) 

On the lattice, the matrix elements of the four-fermion operators above are typ- 
ically extracted from the large-time asymptotic behaviour of three-point correlation 
functions of the form (K p (xi)0 AS=2 (0)K p (x2)} , where K p are pseudoscalar sources 
with suitable quark flavour, K p (x) = d A {x)^s A {x). Expressed in terms of traces 
of quark propagators, these correlation functions correspond to the so-called "eight" - 
shaped quark diagrams given in fig. 

In order to obtain the physical amplitudes, it is necessary to compute the ma- 
trix elements of the renormalized operators corresponding to those defined in eqs. (|2|) 
and @. With Wilson- like fermions, the renormalization procedure is complicated by 
the presence of explicit chiral-symmetry breaking {\SB) in the lattice fermion ac- 
tion: because of the Wilson term, dimension-six operators belonging to different chiral 
representations can mix with each other. 

The testing ground for the restoration of chiral symmetry has been the chiral be- 
haviour of (K°\0 AS=2 \K°) , which, if properly renormalized, vanishes when the K- 
meson becomes massless Q. Although several attempts with Wilson fermions |3|]— |8j 
have given reasonable measurements of Bk, it remains true that the control of the 
renormalization of the relevant operator is rather problematic []. 

1 Bk has also been obtained with staggered fermions Jl^l mainly in the quenched approximation (see 
ref. |]l3f for a review). The (surviving) chiral symmetry in the staggered fermion formalism ensures the 
vanishing of the relevant matrix element in the chiral limit. 
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Figure 1: Connected and disconnected eight diagrams. The dashed line represents the four- 
fermion operator. 



The root of the problem is the operator subtraction outlined above. O^S=2 m ^ xes 
with other operators Oj of the same dimension but with "wrong naive chirality". Thus, 
the (//-dependent) K°-K° matrix element of the renormalized operator 0^ s=2 is given 
in terms of the (a-dependent) bare matrix elements by: 



(i? |O A5=2 (/i)|K c 



lim (K { 

a->0 



I 7 AS=2 
l Z 



AS = 2 (a)+J2M9 2 a)Oi(a) 



\K°) 



The overall renormalization constant Zq S 2 



(4) 

(a/i, g^) is logarithmically divergent and its 
determination does not affect the chiral behaviour of the matrix element. The constants 
Aj^oJ's are finite mixing coefficients which depend on the lattice bare coupling (also 
expressed as (3 = 6/<7g) only [|l4|. They have been either computed in perturbation 
theory (PT) |5|-|lfJ, or fixed by using the non-perturbative (NP) method of refs. p?]]- 
[21 1, or determined using the relevant Ward identities (WI) on quark states (QWI) 



When PT is used for the calculation of the mixing coefficients Aj, the mass de- 
pendence of the renormalized-operator matrix element, (K°\0 AS=2 \K°) , shows large 
deviations from the expected chiral behaviour. Past experience suggests, instead, that 



2 Another method, suggested in refs. [ [23| , 
method has never been implemented though. 



is the use of gauge-invariant WI on hadronic states. This 



the chiral behaviour of (K \O^ s=2 \K°) is satisfactory if a non-perturbative method, 
either the NP renormalization [19]-[2C] or the QWI p2]| , is used. 

In spite of this progress, the determination of the mixing coefficients is a long and 
painful procedure. The non-perturbative renormalization techniques proposed in the 
past have their limitations and, for this reason, an accurate and systematic study of 
the uncertainties in the determination of the mixing coefficients would be necessary 
before keeping the final error on the matrix element fully under control. 

In this paper, we propose a new method which allows the calculation of AS = 2 
amplitudes with Wilson fermions without determining the mixing coefficients. The 
method is based on the lattice Ward identities and it is non-perturbative. It can be 
applied to any Wilson-like formulation of the action (Wilson action, tree-level improved, 
NPI and alike). Its accuracy, in the absence of improvement (of the action and of the 
operators), is ofO(a), which is the best possible accuracy, attainable only with a perfect 
determination of the mixing coefficients Aj in previous approaches. 

The same method can be used to compute AI = 3/2 K-tt matrix elements, e.g. 
those relative to the electro-penguin operators, and AI = 3/2 parity-conserving am- 
plitudes in hyperon decays. It fails, unfortunately, for AI = 1/2 transitions because of 
the mixing with lower-dimension operators. 

An alternative approach, based on the twisted mass QCD formalism pq] , also allows 
the determination of the AS* = 2 mixing amplitude without lattice subtractions pg] . 



2 Description of the Method 

In this section, we show how it is possible to compute the physical K°-K° mixing 
amplitude without knowing, or determining, the mixing coefficients induced by the 
explicit x^B of the lattice action. We first recall some basic notions about operator 
renormalization and mixing with Wilson fermions and then explain our proposal. 



2.1 Operator Renormalization with Wilson Fermions 

Following |2]J (see also || and [27]), we classify the complete basis of dimension-six, 
four-fermion operators which mix under renormalization, relying on general symmetry 
arguments based on the vector-flavour symmetry, which survives on the lattice. To this 
purpose it is convenient to consider separately the parity-even and parity-odd parts 
of the operators in eqs. (Q) and (|3|). Thus, for example, the parity-even part of the 
operator 0\ = Qi — Q\ is given by 

Q x = [VV+AA] = (s A ^d A )(s B ^d B ) + (s A w 5 d A )(s B ^ 5 d B ) , (5) 

whereas the parity odd is defined as 

Qi = [VA+AV] = 2{s A 1 ^d A ){s B 1 ^d B ) . (6) 
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On the basis of CPS symmetries, it can been shown that the renormalization of the Qi 
operators is not affected by the explicit xSB of the lattice action and proceeds exactly 
as in the "continuum" theory. The corresponding renormalization matrix Zij (which 
obviously depends on the renormalization-scheme) is a block diagonal matrix []2l] : 

f Qi\ ( Zn 

Zil Z23 

Z33 





Q 2 
Q 3 
Q 4 



V 











^22 

%32 










Z44 

-2-54 



\ 



^45 
^55 ) 



( Si \ 

0.2 
0.3 

Q 4 

V 0, ) 



(7) 



Thus the lattice does not induce extra subtractions (Aj) for the parity-odd sector (Q&; 
k = 1, . . . , 5), since the mixing in eq. (0) is the same as it would appear in the absence 
of xSB. In particular, Q\ renormalizes multiplicatively. 

The parity-even sector, on the contrary, is not protected by CPS symmetries, and 
all the five relevant operators get mixed because of the xSB of the lattice action. In 
this case, it is convenient to separate the operator mixing into two classes: i) the first 
which consists in correcting the operator mixing induced by the breaking of chiral 
symmetry; ii) the second is the renormalization which survives in the continuum limit. 
In the absence of explicit x^B, the mixing structure is the same as the one considered 
above for the parity-odd counterparts. The corresponding parity-even operators 
k = 1, . . . , 5) would renormalize according to: 

{ Qi\ 



Q 2 
Q3 

V Q5 J 



( Z U 




V 












■^22 


■^23 





■^32 


■^33 











Z44 








Z54 



\ 



Z45 
Z55 J 



( Qi 
Q2 
Q3 
Q4 
V Q 5 



(8) 



where the Qi represent the bare operators, which transform as elements of irreducible 
representations of the chiral group (obviously up to terms of 0(a)). In the presence of 
the Wilson term, the Qi are defined as 

A12 



A 42 
A52 



/ Qi \ 

Q2 
Q3 
Q 4 
V Q 5 / 



/ Qi \ 

Q2 
Q3 
Qa 
Q5 J 



+ 



( 

A21 
A31 
A41 

V a 5 i 



A13 




A43 
A53 



A14 

A 24 
A34 





A15 \ 
A 25 
A35 


J 



I Qi \ 
Q2 
Q3 
Q4 
\ Q5 J 



(9) 



In other words, first the lattice subtraction is performed, followed by the renormaliza- 
tion of the remaining logarithmic divergencies. The above mixing pattern is abbrevi- 
ated, in matrix form, as 

Q = ZQ 

Q = [I + A]Q (10) 

where / is the 5x5 unit matrix. 
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2.2 K°-K° Mixing without Subtractions 

For the sake of illustration, we discuss here the determination of the K°—K° matrix 
element of the operator AS=2 = 0\ only. The extension to the other operators (Oj, 
i = 2, . . . 5) is straightforward. 

Le us consider the Ward identities which can be derived from the T3 axial rotation 

Su = 75U , 5u = ^75 , 
6d = -j 5 d, Sd = -dj 5 , (11) 

where u and d are the up and down quarks taken with degenerate masses, m u = = 
m. For further use, we also introduce the following bilinear operators [] 

U°(x) = d(x)-y 5 d(x) - u(x)~f 5 u(x) , K%{t) = ^ d(x, t)y 5 s(x, t) , 

K° s (t) = J2d(x,t)s(x,t), (12) 

X 

and the corresponding renormalized quantities Kp(t) = ZpKp(t) and Kg(t) = ZgK^it). 
The useful Ward identity in our case is then [|14|] 

(5 [k° p (t^Ox^Kp (i 2 )] ) - (5 [S] Kp^dmkp^)) = , (13) 

where 5 [. . .] denotes the rotation of the argument of 5 and 5 [S] is the rotation of the 
action under the transformation defined in eq. ([ll]) Q 

In terms of the fields defined in eq. (|l2|), and of the parity-even and parity-odd 
operators Q\ and Qi, the Ward identity reads E 



2{k°p(t 1 )Q 1 (0)kp(t 2 )) = 2mJ2(rf(x)KUh)Qi(0)K P (t 2 )) 

X 

-{k° s (ti)Qi(0)k P (t 2 )) - (Kp^Q^kKt,)) + 0(a) , (14) 

where, see eqs. 

Qi = Z n [Qi+ J2 A ^Qi ) . Qi = z nQi- (15) 

The term on the l.h.s. of eq. (|l4|), corresponding to the rotation of the operator Qi, 
is the quantity from which we may extract the physical K°-K° mixing amplitude. 



3 As kaon source, we may as well use the fourth component of the axial current CZ7075S instead of the 
pseudoscalar density. 

4 In general, different choices of the flavour content of the operators, sources and rotation are possible 
or may be necessary. For example for the operators O4 and O5, the suitable choice is O^^ = (sTu)(sTd) {Y 



matrices are not specified) with K\ = d^s and K 2 — U75S) and the same rotation as in eq. (11) 
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Figure 2: Connected and disconnected emission diagrams. The dashed line represents the 
four-fermion operator. We denote as DE{t\,t 2 ) and CE(ti,t 2 ) the diagrams where the K° 
meson is emitted in t\; DE(t 2 , ti) and CE(t 2 , ti) the diagrams where the K° meson is emitted 
in t 2 . 



The first term on the r.h.s. corresponds to the insertion of the rotation of the ac- 
tion. In terms of Feynman diagrams, it can be seen as the "decay" of a neutral (zero 
four- momentum) pion, 7r°, into two K°s under the action of Q\. In the SU(2) isospin 
symmetric case (m u = m^), only the emission diagrams shown in fig. ^ must be con- 
sidered. The last two terms in eq. ( |14| ) correspond to the rotation of the pseudoscalar 
kaon sources. These terms are necessary to saturate the Ward identity. 

One could envisage the following method to extract the coefficients Aij and the 
ratio An = Zn/Zu, thus determining the subtracted operator Q\. Since eq. ( |l~4"l) 
must be satisfied for any value of t\ and t 2 , the Ward identity corresponds to a system 



of linear equations in the unknown quantities Aj (i = 11, 12, . . . , 15) [24]. At least 
in principle, one has an independent equation for any assigned value of t% and t 2 . In 
practice, the equations become dependent when the Ward identity, for large values of 
t\ and t 2 , is saturated by the contribution of the lowest lying states, namely the two 
i^-mesons. 
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This method to extract the Aj, which in practice may be very difficult to imple- 
ment, is however unnecessary and a much easier procedure gives us directly the wanted 
quantity, namely the correlation function of Q\. Let us consider the Ward identity ( |i"4|) 
in the limit t\ —* oo and t 2 — > — oo (in practice for large values of the time distances). 
In this limit, we may safely neglect the last two terms of eq. (|l~4]), because they cor- 
respond to the propagation of scalar states, which are exponentially suppressed with 
respect to the kaon contribution Q Then, up to exponentially suppressed terms, we 
have (An = Z\\jZ\\ and we have divided all the terms of the Ward identity by a 
factor of two) 

t lip A n (K P (h)Qi(0)K P (t2)) -» 

tl— >oo,i2— oo 



\m p\K °)\ 2 



e -m K {t 1+ \t 2 \) x z^{K°\Qi\K°) ~ 



mJ2(n°(x)K P (ti)Qi(0)K P (t 2 )) + 0(a) . (16) 

X 

Note that the last correlation function, (Il (x)K p (ti)Qi(0)K p (t 2 )), is expressed in 
terms of bare quantities only. Indeed a single constant is sufficient to obtain the 
physical amplitude, namely Z\\ which relate the bare lattice parity-odd operator to 
the continuum one, renormalized in a specified renormalization scheme. Obviously this 
constant cannot be determined from the Ward identity, since its values depend on the 



renormalization condition imposed to the renormalized operator |l4| , 21, 22] 



as 



In terms of the Feynman diagrams defined in figs. |T| and [2|, eq. ( |l~6|) can be written 

2{C8(t u t 2 ) + D8(t 1 ,t 2 )) = (17) 
2m(CE(t 1 ,t 2 ) + CE(t 2 ,h) + DE(h, t 2 ) + DE(t 2 ,ti)) + O(a) . 

In summary, the strategy to obtain the physical matrix element (K°\Qi\K°) is 
extremely simple: 

• one computes the correlation function 

G 3 (t h t 2 ) =mY / (n°(x)K P (t 1 )Q 1 (0)K P (t 2 )) 

X 

of the bare parity-odd operator at large time distances ti 2 ; as "mass" m, it is 
more convenient to use the quark mass defined as m = Zap, where p is defined 



using the axial Ward identity [14, ^8], |2^[ 



2p 



(P) 



5 This point can be explicitly checked by computing (K® (t i ) Qi (0) K p (t 2 )) in the same numerical simu- 
lation as the other correlation functions appearing in eq. (|14|). 
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In the above equation and P are the bare (eventually improved) axial current 
and pseudoscalar densities and the matrix elements are usually taken between 
the vacuum and a pion (in our case the tt°) at rest. Za is the axial-current 
renormalization factor. 



then one divides G^(ti,t2) by the factor 

>K 



Ami 



where Z b = \(0\K P \K°)\ 2 , obtaining the quantity K = Z^{K°\Qi\K ). can be 
readily computed from the study of the kaon two-point correlator (K p (t)Kp (0)). 
Alternatively, one may construct the ratio 



\KUh)K°J(0))(K P (t 2 )K P \0)) 



at large time distances. 

finally, the physical amplitude is given by 

{K^Q^K ) = Z n x 11. 

The constant Z\\ cannot be determined from the Ward identity and has to be 
fixed either using perturbation theory, or non-perturbatively on quark states 0]- 
21] or with the Schrodinger functional method [29]. 



The application of the method discussed above to the operators of eq. (jH), appearing 
in extensions of the Standard Model, is so easy that it does not require further discus- 
sion. The same approach can be used to compute the AI = 3/2 K-tt matrix elements, 
{tt\0^ s=1 \K) . A further interesting application is related to hyperon decays. In hy- 
peron decays both the parity-odd and parity-even terms contribute. For the parity-odd 



case, as explained in |2.1| , there are no subtractions induced by the chiral symmetry of 
the lattice action; for the parity-even case, as done before, we add to the relevant 
correlation function a soft-pion (zero four momentum) field, i.e. a pseudoscalar den- 
sity summed over x, multiplied by the factor 2m and replace the parity-even operator 
with the corresponding parity-odd one. Unfortunately, this method works only for 
AI = 3/2 transitions since, in the AI = 1/2 case, the presence of power divergences, 
due to mixing with lower dimensional operators, makes the strategy proposed in this 
paper almost impossible to implement in practice. 

3 Conclusion 

With Wilson-like fermions, the calculation of the K°-K° amplitude is complicated by 
the necessity of accurately determining many mixing coefficients (either perturbatively 
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or non-perturbatively), which arise from the explicit chiral symmetry breaking in the 
lattice action. In this paper, we have shown that, by using suitable Ward identities, it 
is possible to compute the physical amplitudes without any subtraction. The error on 
the final answer is of 0(a), which corresponds to the best accuracy attainable with a 
perfect determination of the mixing coefficients in previous approaches. The extension 
of this method to the improved case, i.e. to obtain matrix elements with an accuracy 
of (D(a 2 ), does not seem possible at present, but it is worth being investigated. Our 
method can also be applied to Al = 3/2 hyperon decays. A pioneering numerical 
calculation of the K°-K° matrix elements using the strategy proposed in this paper is 
already underway. 
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